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Abstract 

The nonperturbative analysis of the one-particle excitation of the electron-positron field is made in 
the paper. The standard form of quantum electrodynamics (QED) is used but the coupling constant 
Qo is supposed to be of a large value (ao S> 1)- It is shown that in this case the quasi-particle 
excitation can be produced together with the non-zero scalar component of the electromagnetic field. 
Self-consistent equations for spatially localized charge distribution coupled with an electromagnetic 
field are derived. Soliton-like solution with a nonzero charge for these equations are calculated 
numerically. The solution proves to be unique if the coupling constant is fixed. It leads to the 
condition of charge quantization if the non-overlapping n-soliton states are considered. It is also 
proved that the dispersion law of the soliton-like excitation is consistent with Lorentz invariance of 
the QED equations. 
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INTRODUCTION 



The problem of quantum electrodynamics (QED) with strong a coupling between 
electron-positron and electromagnetic fields (the so-called "large-a QED") was de- 
veloped by many authors (for example, (l| and citation therein). QED is a physical 
system with well defined Hamiltonian, so the approbation of various methods for this 
theory is of great interest for non-perturbative analysis of actual problems of quantum 
chromodynamics (for example, [i],|3] and citation therein). 

One of the main approaches in this field is related to the Schwinger-Dyson equation 
(SD). the critical value of the coupling constant a c ~ 1 separating the regions of weak 
and strong coupling was discovered for this equation The existence of the critical 
point a c could affect some observed physical phenomena as some authors considered 
p|. Actually the SD equation defines the self-energy part of the fermion propagator 
and results from the accurate QED equations if the corrections for the vertex function 
and the Green's function of the electromagnetic field are not taken into account. This 
corresponds to one of the variants of "quenched QED"when the truncated Fock basis is 
used for the variational perturbation theory A detailed analysis of nonperturbative 
approaches in large-a QED has been recently discussed in paper 

It is important to pay attention to the fact that the mathematical structure of the 
interaction operator in QED is similar in many respects to the operator of the in- 
teraction between electron and quantum field of optical phonons in the ionic crystal, 
the so-called problem of a large-radius "polaron"described by the Frohlich Hamilto- 
nian 0| • This situation was recently used by authors of ji| for nonperturbative mass 
renormalization in the framework of "quenched QED". They made the variational 
estimation of the functional integral for QED on the basis of Feynman's approach to 
the "polaron" problem 

The main effect of the strong interaction between the electron and quantum field is 
defined by formation of the soliton-like ("polaronic") self-localized state of the particle 
coupled with the classical components which are selected from the quantum field by the 
particle itself. It was first discovered by Pekar and investigated rigorously in the 
papers It was also proved that the soliton-like excitation (SLE] corresponds 

to the leading term in the series of a -1 for the Frohlich Hamiltonian [jjj. Later a lot 
of different forms of the quasi-particles arising as a result of a "polaronic "effect were 
described for many concrete physical systems. 

Then the following question is of interest: could some SLE exist in the framework 
of "large-a QED"? It is important to mention that the answer to this question can- 
not be found on the basis of SD equation. The matter is that the integral equation 



for the self-energy function was introduced for the M polaron"problem by Pines 
and it is analogous to SD equation for QED. However, it is well known that 
Pines equation does not allow one to describe the strong coupling limit for the "po- 
laron"problem correctly because the bound state of the electron and phonons is not 
included in this equation. It appeares that the essential reconstruction of the wave 
function both for electron and phonon field should be taken into account in zeroth 
approximation directly if the series in the power of a~ l is used for describing of SLE 
13, 0. If the same problem is considered on the basis of the variational estimation 
for functional integrals some model system qualitatively different from noninteracting 
particles should also be used for calculating the action in zeroth approximation di- 
rectly Particularly, the harmonic interaction between electron and some fictitious 
particle was used for the "polaron"model, that the coupling constant and the mass of 
this particle being considered as the variational parameters 

In this paper it is shown for the first time that SLE can also be found for QED if the 
original coupling constant for the interaction between electromagnetic and electron- 
positron fields is supposed to be large ao ^> 1. As distinct from "polaron"the self- 
localized SLE for QED is the quasi-particle excitation of the electron-positron field 
corresponding to a spatial charge distribution with a non-zero integral charge coupled 
with a classical scalar component of the electromagnetic field separated by the same 
charge density (see also In the terms of "quenched QED"it corresponds to the 

truncated Fock basis which consists of the linear superposition of one-particle exci- 
tation of the electron-positron field, non-zero scalar field and vacuum state for the 
transversal photons of the electromagnetic field. 

The truncation of Fock space used in this paper differs greatly from other non- 
perturbative approaches (see, for example, |l| and references therein) because of the 
inclusion of a strong localized scalar field. However, it should be noticed that a similar 
variant of the quenched QED was recently considered in the paper on the basis of 
mean-field approximation. Some general theorems were proved in this paper but the 
calculation of some concrete characteristics of the system was not given. 

It is also shown in the paper that SLE with a non-zero charge corresponds to the 
minimum of the energy for one-paricle excitation of the system. The characteristics of 
SLE are defined by self-consistent equations for electron-positron wave functions and 
a scalar field. The equations are derived on the basis of an operator method (OM) 
introduced earlier for non-perturbative description of quantum systems ( for example, 
jisl ] and references therein). These results are described in § 2. 

The solution of these nonlinear self-consistent equations is considered in § 3 on the 
basis of the numerical method developed earlier for the "polaron"problem It is 
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shown that a non-trivial solution of these equations exists for a unique value of SLE 
integral charge. It can be considered as the condition of a charge quantization similar 
to that which appeared due to Dirac's monopole ji^l . 

The dependence of SLE energy E(P) on its total momentum P is analyzed in § 3. It 
is well known (2l| that this dispersion law for any one-particle excitation should be of 
a relativistic form because of Lorentz-invariance of QED. The problem of the relation 
between the localization of electron and translational invariance of the whole system 
also appeared for "polaron". It was rigorously solved by authors They proved 

that the "polaron "state vector did not actually break the translational symmetry of the 
system because the concrete point of its localization could be arbitrary in the space [22]. 
As is known, the dispersion law is rather complicated for the "polaron "problem and it 
is very difficult to find its general form analytically because the variables describing 
the internal dynamics of the system are not separated from the translational variables 

The calculation of E(P) for QED (§ 3) showed that SLE localization did not con- 
tradict Lorentz-invariance of the whole system. It is also important that SLE internal 
degrees of freedom are separated from the translational motion and as a result its 
dispersion law E(P) corresponds exactly to conventional dynamics of a relativistic 
free particle. 



II. SELF-CONSISTENT EQUATIONS FOR ONE-PARTICLE EXCITATION IN QED 

Let's consider an exact Hamiltonian QED in the Coulomb gauge j^, This 
representation is most convenient for us because it extracts the electrostatic interaction 
which, by definition, makes the main contribution to the vacuum polarization in the 
strong coupling limit: 



H = I dr{^(r)[a{p+ e A(f)) + (5m^{r) + e o 0(r)p(r) - -(W(r)) 2 } + 
-' 2 

fcA 

Pif) = ^W{r)ijj{f) - %jj{f)^{r)]. (1) 

We suppose here that the field operators are given in the Schrodinger representation, 
the spinor components of the electron-positron operators being defined in the standard 
way |2j] 
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ti(f) = E / (^{4^- e ~^ + V^-a^l- (2) 

In these formulas a, are Dirac matrixes; up sv and ^ Sl/ are the components of the 
bispinors corresponding to the solutions of Dirac equation for the free "bare"electron 
and positron with the momentum p and spin s; a^ s {a^ s ) and bf s (b^ s ) are the annihi- 
lation (creation) operators for the M bare"electrons and positrons in the corresponding 

states. The field operator A(f) and the operator of the photon number n^ A are related 
to the transversal electromagnetic field and their explicit form will be written below. 
It should be remembered that the parameter eo was introduced as a positively defined 
value. 

We should also discuss the dependence of the QED hamiltonian on the operator of 
the scalar field 



<p{f) = V4^ J dkip^. (3) 

It is known [2^|, that these operators can be excluded from the Hamiltonian in the 
Coulomb gauge. For that purpose one should use the solution of the operator equa- 
tions of motion for (p^ assuming that the "bare"electrons are point-like particles and 
"self-action" is equivalent to the substitution of the initial mass for the renormalized 
one. As a result the terms with scalar fields in the Hamiltonian are reduced to the 
Coulomb interaction between the charged particles. However, this transformation of 
the hamiltonian (DO) cannot be used in this paper because just the dynamics of the 
mass renormalization is the subject under investigation. 

There is another problem connected with a negative sign of the term corresponding 
to the self-energy of the scalar field. If the non-relativistic problems were considered 
then the operator of the particle kinetic energy would be positively defined and the 
negative operator with the square-law dependence on (p{f) would lead to the "fall on 
the center" jijj as the energy minimum would be reached at an infinitely large field 
amplitude. However, if the relativistic fermionic field is considered then the operator 
of the free particle energy (the first term in formula (OQ)) is not positively defined. 
Besides, the states of the system with the negative energy are filled. Therefore, the 
stable state of the system corresponds to the energy extremum (!) (the minimum one 
for electron and the maximum one for positron excited states). It can be reached at 
the finite value of the field amplitude (see below). The same reasons enable one to 



successfully use the states with indefinite metric |2a| in QED although it leads to some 
difficulties in the non-relativistic quantum mechanics. 

According to our main assumption about the large value of the initial coupling 
constant eo we are to realize the nonperturbative description of the excited state 
which is the nearest to the vacuum state of the system. The basic method for the 
nonperturbative estimation of the energy is the variational approach with some trial 
state vector |$i > for the approximate description of the one-particle excitation. The 
qualitative properties of the self-consistent excitation in the strong coupling limit IB| 
show that such trial vector should correspond to the general form of the wave packet 
formed by the one-particle excitations of the "bare"electron-positron field. Besides, 
the effect of polarization and the appearance of the electrostatic field (p{r) should be 
taken into account, so we consider |$i > to be the eigenvector for the operator of 
the scalar field. Now, let's introduce the following trial vector depending on the set 
of variational classical functions U^ s ; V^ s ; (p(f) for the approximate description of the 
quasi-particle excited state of the system: 



|$! >~ \<f>?\Uf s ;Vf s ;<p(r)) >= J dq{U^a±+V^b%}\^-ip(T) >; 
£(f)|0;0;<p(f) >=^(r)|0;0;^(r) >; fl? .|0; 0; <p{f) >= 6^0; 0; ^(r) >= 0. (4) 

The ground state of the system is |$o >= |0; 0; >, if we use the same notation. It 
corresponds to the vacuum of both interacting fields. 

Firstly, let's consider the excitation with the zero total momentum. Then the con- 
structed trial vector should satisfy the normalized conditions resulting from the defi- 
nition of the total momentum P and the observed charge e of the "physical"particle: 

< $S 0) |P|$S 0) >= £ dqq\\U<f s \ 2 + |^ s | 2 ] = P = 0; 

s 

E dq\\U qs \ 2 +\V qs \ 2 ] = l. (5) 

s 



< $i 0) |Q|$i 0) >= eoE «*| 2 " l^ s | 2 ] = c (6) 

s 

The condition (0) requires that the functions U qs and V qs should depend on the mod- 
ulus of the vector q only. Besides, one should take into account that the trial vector 

|$i > is not the accurate eigenvector of the exact integrals of motion Q and P as 
it represents the accurate eigenvector of the hamiltonian |$i > only approximately. 
Therefore, in the considered zero approximation the conservation laws for momentum 
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and charge can be satisfied only on average, and this leads to the above written nor- 
malized conditions. Generally, equation (JHj) should not be considered as the additional 
condition for the variational parameters but as the definition of the observed charge 
of the "physical"particle at the given value of the initial charge of the "bare"particle. 
Therefore the sign of the observed charge is not fixed a priori, the same was the case 
for the qualitative estimation in Calculating the sequential approximations to 
the exact state vector |$i > (see §4) should restore the accurate integral of motion as 
well. An analogous problem appears in the "polaron" theory when the momentum con- 
servation law was taken into account for the case of the strong coupling (for example, 
0,0). 

In this respect the variational approach differs greatly from the perturbation theory 
where the zero approximation for one-particle state is described by one of the following 
state vectors: 



|$! >~ \<§>X 1} e >= atjO;0;0 >; |$i >~ \$\ >p >= fctjO; 0; > . (7) 

These vectors don't depend on any parameters and are eigenvectors of the momen- 
tum and charge operators. But they correspond to one-particle excitations determined 
by the charge eo of the "bare"electron and the field cp(r) = 0. We suppose that the 
introduction of the variational parameters into the wave function of the zero approxi- 
mation will enable us to take into account the vacuum polarization, but the variational 
approach in this form allows to conserve the exact integral of motion only on average. 
This problem will be considered in more detail below (§ 4). 

— * 

So, the following variational estimation for the energy E\(P = 0) of the state corre- 
sponding to the "physical"quasi-particle excitation of the whole system is considered 
in the strong coupling zero approximation: 



£i(0) ~ E[ 0) [U qs ; V qs ; <p(f)] =< $<°>| #|$<°> >, (8) 

where the average is calculated with the full hamiltonian (JH) and the functions 
Uqs] Vq S ; ip(r) are to be found as the solutions of variational equations 

dE{°\u qs ;V qs ;cp(r)) = dE? = dE[ 0) = 

dU qs 0V qs dcp(r) [ ' 

under some additional conditions 

It is quite natural, that the ground state energy is calculated in the framework of 
the considered approximation as follows: 
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E ~ 4 0) =< <P \H\<$> > . (10) 

Then the energy of the quasi-particle representing the "physical"electron (positron) 
is calculated by the formula 



E(P = 0) = E{ 0) - 4 0) - (11) 

Further discussion is needed in connection with the application of the variational 
principle ((HIED for estimating the energy of the excited state, because usually the 
variational principle is used for estimating the ground state energy only. As far as we 



know it was first applied in 27], |20| for nonperturbative calculation of the excited 
states of the anharmonic oscillator with an arbitrary value of anharmonicity . This 
approach was called the "principle of the minimal sensitivity". It was shown in our 
paper j^, that the application of the variational principle to the excited states is 
actually the consequence of the fact that the set of eigenvalues for the full Hamiltonian 
doesn't depend on the choice of the representation for eigenfunctions. As a result, the 
operator method (OM) for solving Schrodinger equation was developed as the regular 
procedure for calculating further corrections to the zero-order approximation. Later 
OM was applied to a number of real physical systems and proved to be very effective 
when calculating the energy spectrum in the wide range of the Hamiltonian parameters 
and quantum numbers and the cited references). These results can serve as the 
background for the application of the variational estimation to the excited state energy 
in QED. 

The average value in Eq. (JEj) is calculated neglecting the classical components of 
the vector field (its quantum fluctuations should be taken into account in the further 
approximations) 



< $S 0) |^(r)[aI(f)]^(f)|$S 0) >= 0. (12) 

It should be noted that the possibility of constructing self-consistently the renormal- 
ized QED at the non-zero vacuum value of the scalar field operator was considered 
before ji^l but the solution itself of the corresponding equations was not discussed. 

It is also essential that the normal ordering of the operators wasn't introduced in 
the hamiltonian (JH). Therefore, its average value includes the infinite values corre- 
sponding to the vacuum energies of electron-positron field for both the ground and 
the excited states. However, if the normalized conditions © are fulfilled these values 
are compensated in expression (111ft for the quasi-particle energy. Then the functional 
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for determining the zero approximation for the energy of the one-particle excitation 
is defined as follows: 



E(P = 0) = J drj J0^j2 J ^t) 3 / 2 ^Jl^'s^h'^^ + @ m o)^ + 

e (p(r)5^]U qs u^ - 
Vq's , Vf a , li [{aq + Pmo)^ + e ^(r)^]Vj^}e' (9_9 "' )r - 

\ J driver)} 2 . (13) 

In order to vary the introduced functional let us define the spinor wave functions (not 
operators) which describe the coordinate representation for the electron and positron 
wave packets in the state vector |$[ ^ >: 



In particular, if the trial state vector is chosen in one of the forms (jZj) of the standard 
perturbation theory, the wave functions (THj) coincide with the plane wave solutions 
of the free Dirac equation. For a general case the variation of the functional ( EI ) by 
the scalar field leads to 

f 1 
E(0) = / dr{^ + {r)[{-iaV + (3m ) + -e (p{f)]^{r) - 

y +c (r)[(-iaV + (3m ) + \e Q <p(f)]^ c {r)- 

J dr[^ + (f)^(r) + ty +c (r')ty c (f')} = 1; (15) 



^(r) = ji/ p^[* + (r~W) " ^ +c (^)^)]; (16) 

| df[*+(r)*(r) + ^ +c (f )^ c (r f )} = 1. (17) 

The main condition for the existence of the considered nonperturbative excitation 
in QED is defined by the extremum of the functional ( fT5| ) corresponding to a non- 
zero classical field. The structure of this functional shows that such solutions of the 
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variational equations could appear only if the trial state vector simultaneously included 
the superposition of the electron and positron wave packets. 

Equation f fT5| ) and the Fourier representation (J3J) for the trial vector clearly indicate 
that the assumption concerning the localization of the functions \I/(r) near some point 
doesn't contradict the translational symmetry of the system because this point can 
itself be situated at any point of the full space with equal probability. The general 
analysis of the correlation between the local violation of the symmetry and the con- 
servation of accurate integral of motion for the arbitrary quantum system was first 
discussed in detail by Bogoluibov in his widely known paper "On quasi-averages"|22]. 
A similar analysis of the problem in question will be given in § 4. 

Varying the functional ( TT51 ) by the wave functions \I/(r) h \I/ C (r) taking into account 
their normalization leads to the equivalent Dirac equations describing the electron 
(positron) motion in the field of potential <-p(f): 



{{-iaV + Pm ) + e (p(r)}^(r) = 0; 

{(-iaV + (3m ) + eo(p(r)}^ c (r) = 0. (18) 

But it is important that in spite of the normalized condition (ITT}) for the total state 
vector ( IHl ) each of its components could be normalized differently 



r 1 

J df^ + (r)^(r) = j-^; 

df^ +c (r f )^ c (r f ) = - (19) 

1 + c 

The constant C is arbitrary up to now. It defines the ratio of two charge states in 
the considered wave packet. As a result the self-consistent potential (p{r) of the scalar 
field depends on C because of the equation f TTHD . 

We should discuss the procedure of separating variables in more detail, because of 
the non-linearity of the obtained system of equations for the wave functions and the 
self-consistent potential. Since the considered physical system has no preferred vectors 
if P = 0, it is natural to regard the self-consistent potential as spherically symmetrical. 
Then the variable separation for the Dirac equation is realized on the basis of the well 
known spherical spinors (25| : 



jlM 



g(r)ttjiM 
if(r)VL jVM 



(20) 
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Here QjiM are the spherical spinors describing the spin and angular dependence 
of the one-particle excitation wave functions; j, M are the total excitation momentum 
and its projection respectively, the orbital momentum eigenvalues are connected by the 
correlation l+V = 2j. It is natural to consider the state with the minimal energy as the 
most symmetrical one, corresponding to the values j = 1/2; M = ±1/2; I = 0;l' = 1. 
This choice corresponds to the condition according to which in the non-relativistic limit 
the "large "component of the spinor \I/ ~ g corresponds to the EPF electronic zone. 
Then the unknown radial functions /, g satisfy the following system of the equations: 



- \rg) -(E + mo- e p(r))(r/) = 0; 
ar r 

d ^ + \rf) + (E — m - e ip(r))(rg) = 0. (21) 
ar r 

The states with various projections of the total momentum should be equally popu- 
lated in order to be consistent with the assumption of the potential spherical symmetry 
with the equation ifTBjl . So, the total wave function of the "electronic"component of 
the EPF quasi-particle excitation is chosen in the following form: 

* = ^=[^1/2,0,1/2 + ^1/2,0,-1/2] 
1 



9\r)x6 
if{r)xt 



xt = -^Pi/2,/,1/2 + ^1/2,1,-1/2}; l = 0; 1. (22) 
In its turn, the wave function ^ c is defined on the basis of the following bispinor: 



fC I -ifi(r)ttjiM 
i m \ gi{r)^ Vh 

The radial wave functions /1, g\ in this case satisfy the following system of equations 



\ gi(r)VL jVM 1 



d ^ + -( r 4i) " (E + m + e p(r))(r/i) = 0; 
ar r 

d ^p± - -(rfi) + (E - m + e ip(r))(r gi ) = 0. (24) 
ar r 

These equations correspond to the positronic branch of the EPF spectrum with the 
"large "component ~ g\ in the non-relativistic limit. 

It is important to note that the functions \I/ and \I/ C satisfy the equations f fTKl ) with 
the same energy E(0). It imposes an additional condition of the orthogonality on them: 
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< * c |* >= 0. (25) 

Taking into account this condition and also the requirement that the states with 
different values of M should be equally populated one finds the "positronic"wave 
function 



* C= -7^[*lA0,l/2-*l/2,0,-l/2] 

XT = -7^1/2,1,1/2 - ^1/2,1,-1/2]; I = 0; 1. (26) 



1 r JC JC ] — ( ~ ^/l( r )X0 

^L W l/2,0,l/2 - W l/2,0,-l/2J - I ^( r ) x - 
1 

The equation for the self-consistent potential follows from the definition of (fi(r) in 
formula (fT%I)taking into account the normalization of the spherical spinors (25| : 

d 2( P , 2 di P e o r,2 , 2 ,2 21 

The boundary condition for the potential is equivalent to the normalization condition 
(ED and defines the charge of the "physical"electron (positron) e 

f rldn[f\n) + ff 2 (n) - /,» - gl(n)]. (28) 

It is important to stress that the form of the functions given above is defined prac- 
tically uniquely by the imposed conditions. At the same time the obtained equations 
are consistent with the symmetries defined by the physical properties of the system. 
The first symmetry is quite evident and relates to the fact that the excitation energy 
doesn't depend on the choice of the quantization axis of the total angular momentum. 

Moreover, these equations satisfy the condition of the charge symmetry j^j. Indeed, 
by direct substitution, one can check that one more pair of bispinors leads to the 
equations completely coinciding with (Ell), (Ell) 

= f i9 f^' M ) ; (29) 
\-fi(r)n jlM ) 

**» " { -ig(r)n jlM I ' < 30) 
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It means that these bispinors allow one to find another pair of the wave functions 
which are orthogonal to each other and to the functions (j22j),(J2E|) but include the same 
set of the radial functions 



*=U$)- (31) 

These functions differ from the set (j22J),(J2Ej) in that they have a different sign of the 
observed charge of the quasi-particle due to the boundary condition ( ESI ) - 

Let us now proceed to the solution of the variational equations. It follows from the 
qualitative analysis that the important property of the trial state vector freedom is the 
relative contribution of the electronic and positronic components of the wave function. 
Therefore let us introduce the variational parameter C in the following way: 



/*00 1 

J q r 2 dr[f 2 (r)+g 2 {r)] = — - 



fr 2 dr{fi(r)+g 2 (r)] = T ^. (32) 
The dimensionless variables and new functions can be introduced 



x = rm ; E = em ; e (p(r) = m o 0(x); — = a ; 

u(x)y/m^ = rg(r); v{x)^/rn^ = r/(r); ui{x)^Jrn^ = r#i(r); ui(a;) x /m^ = r/i(r(B3) 

As a result the system of equations for describing the radial wave functions of the 
EPF one-particle excitation and the self-consistent potential of the vacuum polariza- 
tion can be obtained. 



u — (e + 1 — (f)(x))v = 0; 

CL00 00 

dv 1 . „ , , . . 

— + -v + (e - 1 - (f){x))u = 0; 

CL00 00 

du\ 1 

— h —u\ — (e + 1 + <p{x))v\ = 0; 

CLOC 00 

dv\ 1 

— v\ + (e - 1 + <p[x))u\ = 0; 

CLOC 00 

poo \ U) J r% 

<t>{x) = a [ dy h - / dyp(y)}; 

Jx y x JU 

p{x) = [u 2 (x) + v 2 (x) — u\{x) — v\{x)]. 



(34) 
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III. NUMERICAL SOLUTION OF THE SELF-CONSISTENT EQUATIONS 



The mathematical structure of equations Q31D is analogous to that of the self- 
consistent equations for localized state of "polaron"in the strong coupling limit 
(ll|,(l2|. Therefore the same approach can be used for the numerical solution of these 
nonlinear equations. It has been developed and applied |3| for the "polaron" problem 
on the basis of the continuous analog of Newton's method (CANM) • 

Let us take into account that the system of equations (Ell) can be simplified because 
the pairs of the functions it, v and ui, v\ are satisfied by the same equations and differ 
only by the normalized condition. Therefore they can be represented by means of one 
pair of functions if special notations are used: 



u 



\ 1 + C 



uo; v 



1 

\TTc 



C 

\l + C 



vq; vi = 



C 

\TTc U0] 



J o dx[u 2 (x) + Vq(x)] = 1; po(x) = u\(x) + Vq(x); 



dx 
dvQ 



1 — C 

<f)(x) = a - - c <J)q(x); (j) (x) = 



1 

uo 

X 

fir + x V ° ~ ^ + = 0; 

dy^l + - [ dyp (y)}. 
y x J o 



The energy of the system (jTHjl can also be calculated by these functions 

E(o) = m „i T g[r + ia i T §n]; 

T = [ dx[(u v - v' u ) - + (uq - Vq)]; 



n 



X 



dX(f) (ul + Vq). 



(35) 



(36) 



and equations can be obtained as the result of the variation of the functional 



The required solutions are to be normalized and this condition defines the asymptotic 
behavior of the functions near the integration interval boundaries: 



uq pa Ax[l + 



1 



6 



Vq 



A 



1 



X 



0; 



The quantity 



uq « A\ exp(— x), Vq « — A\ exp(— x) x — > oo. 
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(37) 



1 — c 

a = a °Y~- Tc' ^ 

is a free parameter and it should be chosen from the condition of the existence of 
non-zero and normalized solution of equations (ESI)- 

In order to use the method developed in ref.|ll| especially for nonlinear eigenvalue 
problems let us consider instead of (ESI) the following system of equations: 

du k 1 , , 2 



doc oc 



Uk- (k - (f)(x))vk = 0; 



dvk 1 /j2 if w 

~dx + x Vk ~ ' ( f > ( x )) u o = °5 

| Q dx[u\{x) + ^(x)] = 1; p k (x) = u\{x) + 
(f)(x) = a(f) k (x); (f) k (x) = [f dy^^- + - f dyp k (y)}. (39) 

Here the eigenvalue k 2 (a) is introduced when the parameter a is fixed. In this case 
the normalization and the boundary conditions 

, ri A: 2 - 0l(o) A k 2 - (f)(0) 2 

6 3 

Wfc « kAi exp(— kx), Vk ~ — exp(— /ex) x — ► oo (40) 



allow one to calculate k 2 (a) by means of the numerical method from ref. 19]. Then 
one should change the parameter a in order to find the value ao which satisfies the 
condition 



k 2 (ao) = 1- (41) 

Let us consider briefly the iteration procedure for calculating the functions 
Uk{x),Vk{x) and the value k. If (s) iterations are made, the desired values are de- 
fined by some approximate expression 



4 s) (x), v[ s \x) *«; 
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roo 



£ P { :\x)dx = i, p { k s \x) = [i4 s) wi 2 + i4 s) w 

4 S \x) = a[jy y ^± + 1 £dyp K %)\. (42) 



2/ z 

Then the following recurrent equations describe the wave functions and eigenvalues 
for the subsequent iterations: 



u[ s+1 \x) = A^{ui s \x) + r[^ s \x) + ^\x)}}, 
v[ s+1 \x) = A^{v[ s \x) + r[4 s \x) + P is) 4%)}}, 

k ( S +l) = k (s) + ^ (43) 

and the exact solutions are calculated as the limits of sequences 



u k (x) = lim 4 S, W, v k (x) = s lim4 s) W, k = lim /c (s) . (44) 
Here r < 1 is an arbitrary parameter determining the step of discretization for the 



evolutional equation corresponding to CANM |20| for the considered problem 
The corrections i[)( s \x), ip[ (x),i[)^\x), ipi^(x) are supposed to be small if the zeroth 
approximation is good. They can be calculated as the solutions of the nonhomogeneous 
equations resulting from the initial problem jiol ] 

d .fl _ >> _ (fc <* _ ^>)# = _ I U M _ (fc (* - 

CLOO 00 CLOC 00 

CLOC 00 CLOC 00 



(45) 



doc oo 

The system of differential equations should be supplemented by the boundary con 
ditions resulting from ( IHTj ) 

ip\ >(x) fc v j xi/j [s) (x) = -[v {s) (x) fc v J xu (s) (x)}, X -> 0, 

o o 

/c (s Vi S) (^) + ^ (s) 0) = -[£ (s V s) (:c) + u {s) (x)}, x -»■ oo, 
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\ x ) + ^O) = -v {s) (x), x -»■ oo.(46) 

The correction yi/.^ for the eigenvalue /c^ and the normalization amplitude are 
calculated on the basis of the normalization condition for the wave function. The first 
order correction for the normalization integral tends to zero by means of the available 
choice of the value // s ) : 



/„« = / o °°K'Vi 5) + ^Vg]^- (47) 

The coefficient ^4^^ allows one to normalize functions (IIHI) for the subsequent iteration 
with the second order accuracy on the value /i^. The convergence of the considered 
iteration scheme was proved in ref.|3| . 

It is known j^H, that the numerical solution of the differential equations ( 1431 ) within 
the interval [0, oo] is more effective if the following logarithmic change of variables is 
used 



x = e 9 , 9 = lux, dx = e d d9, -oo < < oo, (48) 

It allows one to distribute the nodes for the finite-difference approximation of equations 
( l4*5l ) in a nonuniform way. As a result the step of integration is small in the range 
x < 1, where the desired functions are most significant, and increases for large values 
of x when the functions are small. 

Zeroth approximation for the recurrent equations ( 1151 ) can be calculated by means 
of the variational estimation for the functional (It36j) using the following trial functions 





4 0) = Ax(l + bx)e~ bx , 4 0) = Bx 2 e- hx 1 

7 A 2 + 3B 2 = 46 3 . (49) 

Thus, for example, the recurrent sequence (jUJ) converges quite quickly into the value 
k(a) = 1.05 for the following initial values of the parameters 



A = J2j7, B = ^3, k {0) = 1, r = 0.5, 6=1, a = -3.3. (50) 
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It is important that the variation of the initial approximation for the functions and 
parameters affects the speed of the iteration convergence but does not change the final 
result for k(a). It corresponds to the only global extremum for the functional (jHED- 

The numerical solutions Uk(x), Vk(x), k(a) calculated with the fixed value a were 
used further as the initial approximation for the iteration solution of the equations 
( BSj ) w ^ n t ne variation of the parameter a\ = a + Aa, This procedure was repeated 
until the condition (J2J) was satisfied with some given accuracy. As a result ao was 
calculated up to 4 significant figures: 



a ~ -3.296, (51) 

and this value corresponds to the desired soliton-like solution with nonzero charge. 

The numerical analysis of the equations in wide range of parameters has shown 
that the only solution described above existed under the condition that the function 
wo should not have zero points for the ground state. The value ao can also be expressed 
in terms of the integrals from the functional provided of its minimum relative to 
the parameter 5 = (1 — C)/(l + C), which defines the contributions of the electron 
and positron components to SLE : 

T / /-v \ mnTa . . 

<*o = - W ; E(0) = ~- 52 

11 ao 2 

The calculated value of the integral 



T = 0.749. (53) 

FigjT] shows the solutions uq,vq of the Dirac equation and FigEl represents the self- 
consistent potential 0o corresponding to the value ao- The characteristic size Ar of 
the spatial localization of the excitation is defined by the parameter 



Ar^niQ 1 . (54) 

The equations ( ESI ) and (EH) lead to the following relation between the "bare "charge 
eo and the integral charge of SLE: 



eo 



= e; 



e eg = 47ra, e = 



1 + C) 
47ra 41.42 



eo 



eo 



(55) 
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Phc. 1: Solutions of the Dirac equation for SLE 

It is of interest that the formula (ESlhas the same structure and consequences that 
result from the Dirac monopole theory 20| despite of a different physical interpretation. 
In particular, if the calculated self-localized one-particle excitation of the electron- 
positron field (EPF) could be considered as a "physical "electron the only possible 
value of the parameter ao together with the ratio ( JH51 ) leads to the condition of the 
observed charge quantization with the fixed value of the "bare"charge eo. Indeed, on 
the analogy with the "polaron"theory (s2|,(s3| n-particle excitation of EPF can be 
considered as the superposition of n different SLE' that are at the distances exceeding 
the size of localization Ar. Then the charge of such n-soliton excitation q n would be 
always multiple to the elementary charge of SLE q n = ne. 
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Phc. 2: Self-consistent potential of SLE 
IV. EPF QUASI-PARTICLE EXCITATION WITH AN ARBITRARY MOMENTUM 

In the previous section we have considered the possibility for a resting quasi-particle 
with a non-trivial self-consistent charge distribution, the finite energy E(0) and a zero 
total momentum P = to exist in the framework of a nonperturbative QED. 

The obtained solution allows one to imagine the internal structure of the resting 
"physical"electron (positron) as a strongly coupled state of charge distributions of the 
opposite sign, the large values of integral charges of these distributions compensating 
each other almost completely and their heavy masses being "absorbed "by the binding 
energy. 
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Actually the energy ±£"(0) defines the boundaries of the renormalized electron and 
positron zones resulting from the strong polarization of EPF when the excitation 
appears. But this excitation could be interpreted as the "physical"electron (positron) 
if the sequence of the levels in every zone determined by the vector P / (Fig.l) 
were described by the relativistic energy spectrum of real particles, that is 



E(P) = /P 2 + £ 2 (0) = VP 2 + m 2 . (56) 

Only in this case the energy E(0) can really be used for calculating the mass mo of 
the "bare"electron. 

It is worth saying, that the problem of studying the dynamics of the self-localized 
excitation should be solved for any system with a strong interaction between quantum 
fields in order to calculate its effective mass. For example, a similar problem for Pekar 
"polaron"|Io| in the ionic crystal was considered in papers (H], |3| and 

in many more recent works. It is essential that because of the non-linear coupling 
between the particle and a self-consistent field the energy dispersion E(P) for the 
quasi-particle proves to be very complicated. As the result, its dynamics in the crystal 
is similar to the motion of the point "physical"particle only at a small enough total 
momentum. 

However, in the case of QED the problem is formulated in a fundamentally different 
way. There is currently no doubt that the dynamics of the "physical"excitation should 
be described by the formula (EED for any (!) values of the momentum P. It means that 
the considered nonperturbative approach for describing the internal structure of the 
"physical "electron should lead to the energy dispersion law ( 151)1 ) for the entire range 
of the momentum P. 

The rigorous method of taking into account the translational symmetry in the strong 
coupling theory for the "polaron" problem was elaborated in the works of Bogolubov 
and Gross Let's remember that this method was based on the introduction of 

the collective variable R conjugated to the total momentum operator P, the canonical 
character of the transformation caused by three new variables Ri being provided by 
the same number of additional conditions imposed on the other variables of the system. 
In the "polaron"problem the quantum field interacting with the particle contributes 
to the total momentum of the system. It allows one to impose these conditions on the 
canonical field variables JJ] and the concrete form of the variable transformation 



is based mainly on the permutation relations for the boson field operators. 

The considered problem has some specific features in comparison with the "po- 
laron"problem. Firstly, the formation of the one-particle wave packet is the multi- 
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particle effect because this packet includes all initial states of EPF as the fermionic 
field. Secondly, its self-localization is provided by the polarization potential of the 
scalar field that doesn't contribute to the total momentum of the system. Therefore, 
we use a different approach in order to select the collective coordinate R. Let us return 
to the configuration representation in the Hamiltonian (JTJ) , where QED is considered 
to be the totality of N (N — > oo) point electrons interacting with the quantum EMF 
in the Coulomb gauge [23]: 



N 

H = 2 {a a [p a + eoA(f a )} + /3 a m + e ^(r a )} 

a=l 



df(V^(f)) 2 + E^(/c)% A ; 



fcA 

u(k) = k; % A = ctc^; A = 1,2; 



m = g sm^ x) h^ ,kf + cy- ,tf ] (57) 



Here Q is the normalized volume; ci" (c^ A ) are the operators of the creation (annihi- 
lation) of quanta of a transversal electromagnetic field, the quantum having the wave 
vector k, polarization and energy tu(k) = k. The sign of the interaction operators 
differs from the standard one because the parameter eo is introduced as a positive 
quantity. 

As it was stated above the zero approximation of nonperturbative QED is defined 
only by a strong interaction of electrons with the scalar field, where the interaction 
with the transversal field is to be taken into account in the framework of the stan- 
dard perturbation theory when the conservation of the total momentum is provided 
automatically ji^. Therefore, while describing the quasi-particle excitation we should 
consider the conservation of the total momentum only for the system of electrons. In 
the considered representation it is defined by the sum of the momentum operators of 
individual particles 

N „ 
a=l 

P|$i(P) >= F|$i(P) > . (58) 

It means that in the configuration space the variable R conjugated to the total 
momentum is simply a coordinate of the center of mass and the desired transformation 
to new variables is as follows: 
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1 N N 

r a = R + p a ; R = ^J2r a ; E Aj = °i 

iV a=l a=l 

1 -> ~ % - -> % N -> N a 

Pa = -iV a = --V R + p a , P=-iV R ; f a = -iV Pa + - E V Pb ; E# a = (B9) 

^ v ^ v 6=1 a=l 

The Hamiltonian ( 1571 ) with new variables is of the following form 

N i - 

# = E{«a[-T7Vi? + e A(p a + iQ] + /? a m + e o 0(p a + #)} - 

a=l iV 

i/^r(V^(r)) 2 + E^(^ A - (60) 

It should be noted that the matrix elements of an arbitrary operator in a new con- 
figuration representation are to be calculated in accordance with the following norm 
we introduce a special notation for this norm) 

« $!|M|$ 2 »= / dRUdpa^KR, {p a })M<S> 2 (R, {p a }). (61) 



Let us denote by Ho that part of the operator f lHOl ) which doesn't depend on the 
transversal EMF and describes the internal structure of "physical"particles. In fact, the 
operator Hq doesn't depend on the coordinate R either, because of its commutativity 
with the operator of the total momentum of the system of electrons. This also follows 
from the well known result 0| that in the Coulomb gauge the scalar potential could be 
excluded from the hamiltonian. As a result the operator H depends only on the vector 
differences (f a — r^) = (p a — p&) and doesn't change with the simultaneous translation 
of all the coordinates. As a consequence, the eigenfunctions of the hamiltonian Ho 
depend on the coordinate R in the same way as for a free particle: 

HR: {Pa}) = ^372 e ^l $ l(^ te» >5 

H - H (P) = Eft4p+a + P a m + e o 0(p a )} - \ I dr(V(p(r)) 2 (62) 

Further calculations consist in returning to the field representation by the variables 
p a in the limit (N — > oo) and in using the approximate trial wave packet |$i(P) > 
similar to (jlj) but with the coefficient functions depending on P . Thus, in the frame- 
work of non-perturbation QED the orthogonalized and normalized set of states for the 
EPF one-particle excitation is defined as follows 
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l*S° \P) >^ (^p e ^/ dRUdPK + V&(P)bt,}\0; 0; tp(f) >; (63) 

with the norm (1ETT) and the coefficient functions U^ S (P)] V^ S (P). 
Using the coordinate representation for these functions 



^) = / 7^72 E ^ s (P)^e^ T , (64) 



one can find the following functional for calculating the value E(P) = e[ ^ — Eq 
corresponding to the energy of one-particle excitation with an arbitrary momentum 

E(P) = [ df{^ + (f, P)[(aP - laV + /3m ) + e ~(p(r, P)]*(r, P) - 
^ +c (r, P)[(-aP - iaV + (3mo) + e ^(r, P)]* c (r, P); 

<p(f, P) = e / ^^[^ + (r , P)*^, P) - ^ +c (^, P)*V, P)]; 

| dr[* + (r, P)#(r, P) + * +c (^, P)* c (r", P)] = 1. (65) 

It should be noted that the mean value of the operator Hq(P) in the representation 
of the second quantization was calculated by taking into account the expression 

lim — V = 1, 

1 * p 

that follows from a well known formula for the density of states when EPF is quantized 
within the normalized volume Q (j^ . 

It was mentioned above, that there is quite a close analogy between the nonpertur- 
bative description of QED and the theory of strong coupling in the "polaron"problem. 
Therefore, it is of interest to compare the obtained functional (E5j)with the results 
of various methods of including the translational motion in the "polaron"problem. 
The simplest one was used by Landau and Pekar j^ij who introduced the Lagrange 
multipliers in the form 



j(p) = j(p = o) + (pv), 



(66) 
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with the functional J(P = 0) referring to a resting "polaron"and the Lagrange 
multiplier V{ denoting the components of the quasi-particle average velocity. 

We can see that the obtained functional (E5D has the same form as flfifitlif the rela- 
tivistic velocity of the excitation is determined by the formula 

V = |rff{^+(f,P)(a)^(r,P) + ^ +c (f,P)(a)^ c (f,P)}. (67) 

This corresponds to the well known interpretation of Dirac matrixes. 
Varying this functional by taking into account the normalized conditions leads to 
the following equations for the wave functions 

{(aP - laV + /3m ) + e (p(f, P) - E(P)}V(r, P) = 0; 

{(aP + iaV - (3m ) - e (p(r, P) - E(P)}^ c (r, P) = 0. (68) 

These equations show that unlike the "polaron"problem the translational mo- 
tion of the quasi-particle determining the momentum P is related in our case to its 
internal movement described by the coordinate r by means of spinor variables only. 
The physical reason for this separation of variables is explained by the fact that in 
QED the self-localized state is formed by the scalar field, and its interaction with 
the particle doesn't involve the momentum exchange. In order to find the analytical 
energy spectrum E(P) the system of non-linear equations (JHB1 should be diagonalized 
relative to the spinor variables. The possibility of such diagonalization seems to be a 
non-trivial requirement for the nonperturbative QED under consideration. 

The solution of the equations (jTTHT) can be found on the basis of the states for which 
the dependence on the vector (fin the wave packet amplitudes U$ )S ,V% S remains the 
same as it was in the motionless "physical"electron. However, the relation between 
the spinor components of these functions can be changed. But as the self-consistent 
scalar potential involves the summation over all spinor components it can be assumed 
that the potential doesn't depend on the momentum for the class of states in question: 

<p(r,P) = ip(r)\ p=Q . (69) 

In the coordinate representation the transformation of the spinor components of the 
wave functions satisfying the equations ( 1681 ) takes place because of the dependence on 
the momentum. It seems that there are some general arguments based on the Hamil- 
tonian symmetry that make the diagonalization of equations (jSBj) easier. However, 
we have desired the solution by sorting out various linear combinations of the wave 
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functions found in Section 3 for a resting electron. These functions correspond to the 
degenerated states in the case of P = but are mixed for a moving electron. It is 
found that there is only one normalized linear combination satisfying all the necessary 
conditions of self-consistency : 

#(r, P) = L{P)V{r) + K{P)&{r); 
^ c (r,P) = Li(P)^ c (r) + i^i(P)§(r); 

|L| 2 + \K\ 2 = \L X \ 2 + \K X \ 2 = 1. (70) 

The condition (EDJ) according to which the potential doesn't depend on the momen- 
tum for the excitation at the energy E(P), is fulfilled if the coefficients are related 

as 



Li = -K; K x = L. (71) 

These relations are also consistent with Equations ( IHHl ) for wave functions. 

This means that the "physical"electron moves in such a way that its states are 
transformed in the phase space of the orthogonal wave functions (J22I) , (EEI) , (EDI) , (Ell) but 
the amplitudes of its "internal "charge distributions are not changed. These results, 
however, are valid only for the case of neglecting the interaction with the transverse 
electromagnetic field. 

Substituting the superpositions ( 1701 ) into equations (EBD we use the following rela- 
tions: 



\ g(r)(0-P)it l/2A M J \-!/l('-)(CT-P)fi 1/ 2,o,M 

H5V + Pm + e oV )m = E(0) ( fj;^ M 

\lj[r)\Li/2,l,M 

(_. i5 V + 0m o + e oV W(r) = -E(0) ( ) ,(72) 

and similar formulas for the functions ^(r); \I/ C (r); Oi are the Pauli matrixes. 

For equations ( IHHj l to be fulfilled for any vector Fit is necessary to set the coefficients 
of spherical spinors equal to the same indexes /. The corresponding radial functions are 
proved to be the same under these conditions, and the following system of equations 
for the spinors xti i s obtained 
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iL(aP)xt + K(E + E Q )xt = 0; iK{aP)xt + L(E - £ )x + = 0; 

L(aP)xi + + ^o)Xo + = 0; K(aP)xt + - E )xt = 0; 
L l (aP)xi ~ iKi(E - E ) X i = 0; K l {dP) X , - iLx{E + E ) Xo = 0; 
iL 1 {&P)x[ -ME- E )xi = 0; iK^P) X ^ ~ HE + £ )Xo = 0. (73) 

Spin variables in Equations (1731 are also separated. In order to show this one can 
use, for example, the relation between the coefficients resulting from the 4th equation 
in ( 1731 in the first one of these equations: 

+ = . K(aP)xt 

Xl l L(E-E y 

As aresult there exists a non-trivial solution of these equations for two branches of 
the energy spectrum 



E ejP = ±y/E% + P\ (74) 

referring to the electron and positron zones, respectively (Fig.l). The same expres- 
sions can be obtained for all conjugated pairs of the equations in ( 1731 . The coefficients 
in the wave functions ( 17(1 can be found taking into consideration the normalization 
condition: 



L e = K{= . P ■ K e = -L{ E< E{ 



P 2 + (E e -E )^ ^P* + (E e -E ] 

P r n E g + Eq 
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L p = K{= ; K p = —L\ = — . (75) 

lpi + (E e + E^ JP* + (Ee + E,f y ) 



One can see that this set of coefficients coincides with the set of spinor components 
for solving the Dirac equation for a free electron with the observed mass m = Eq. 
Thus, the results of this section show that the "internal"structure of the "physi- 
cal" electron (positron) considered in this paper is consistent with the experimental 
energy dispersion ( 1741 for a real free particle due to the relativistic invariance of the 
Dirac equation. 

For the interpretation of the wave packet ( 1H31 as the state vector for a "physi- 
cal"particle it is essential that it is the eigenvector of the total momentum of EPF in 
the framework of the considered "quenched QED". It means that SLE with different 
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momenta form the set of orthogonal and normalized functions if the condition ( I61j ) is 
taken into account: 



Thus, in the present paper new variational self-consistent equations for interacting 
electron-positron and scalar electromagnetic fields have been derived in the framework 
of the "large-a QED". These equations differ from the Schwinger-Dyson equations and 
describe the one-particle excitation of the system that doesn't contain the quanta of 
the transversal electromagnetic field. 

The soliton-like solution of the desired equations has been calculated numerically. 
It includes a non-zero integral charge ("charge"soliton) and can be regarded as the 
spatially localized collective excitation of the electron-positron field coupled with the 
localized scalar field. It is shown that this solution is unique if the coupling constant 
ao for the field interaction is fixed. This result leads to the quantization condition for 
the "observed" charge of n-soliton excitations of the system. 

The dependence of the energy of the soliton-like excitation on its total momentum 
has been also studied. It is shown that the corresponding law of the energy dispersion 
is consistent with Lorentz invariance of QED and the effective mass of the "charge 
soliton"can be calculated. 

It should be emphasized that the main aim of the paper is to prove the possibility of 
the existence of such non-trivial solutions of the "large-a QED"equations that could 
not be derived on the basis of the standard perturbation theory. However, the physical 
interpretation of the desired solution requires further research. 
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